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■ The mean field theory of a spin glass with a specific form of nearest and next nearest neighbor 

interactions is investigated. Depending on the sign of the interaction matrix chosen we either 
, find the continuous replica symmetry breaking seen in the Sherrington-Kirkpartick model or a one 

/"*\ ■ step solution similar to that found in structural glasses. Our results are confirmed by numerical 

simulations and the link between the type of spin-glass behavior and the density of eigenvalues of 
r*"** . the interaction matrix is discussed. 
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I. INTRODUCTION 

The Shemngton-Khjkpatrick (SK) modelEI is the most well-known example of a disordered and frustrated system in 
the field of spin glassesa. It corresponds to the infinite-range version of the Edwards- Anderson (EA) model introduced 
earlier in 1975a. In the EA model quenched disorder is introduced in the random sign of the exchange couplings between 
nearest-neighbor spins on a lattice. The infinite-ranged version is the natural mean-field version of the EA model, 
in the same sense as the infinite-range ferromagnet is the mean-field theory for the Ising model. A nearly complete 
solution of the SK model has been founda which has raised subtle questions about the nature of the spin-glass phase 
showing that the mean field theory of spin glasses is considerably more complex than the standard mean field of 
say ferromagnetic systems. Plenty of new questions were posed after it was shown that the correct thermodynamic 
solution had to be understood in terms of a replica-symmetric broken solution. For instance, does the ergodicity 
breaking (implied by replica symmetry breaking) also occur in short-range systems? Also, has replica symmetry a 
I ' true physical meaning in as much as time-reversal symmetry has for usual ferromagnets? After many years of research 
, these questions have turned out to be extremely difficult and in the meantime criticisms questioning their relevance 
\Q • regarding our understanding of experimental systems have also been raised. Scaling theories of finite dimensional 
spin glasses, so called droplet modelsu, seem to be at variance with the image of replica symmetry breaking. It is 
■ thus useful to introduce new solvable models which correspond to mean field versions of different finite dimensional 
' problems in order to improve our understanding of the spin-glass problem. 

In this paper we introduce a new solvable spin-glass model which corresponds to the mean-field version of the 
Edwards- Anderson model but including next-nearest neighbor interactions. We will refer to this new model as NNN- 
SK model. The motivation is that the model admits the possibility of being realized on a finite-dimensional lattice 
and incorporates correlations between the first nearest-neighbor and second- nearest neighbor coupling interactions. 
At a first glance it seems strange to consider a next nearest neighbor interaction in a totally connected spin-glass 
1 ■ model, however we shall see that the correlations introduced in the couplings can lead to a physics different to that 
of the SK model. In addition this model allows one to discuss the role of the density eigenvalues of the interaction 
matrix in the spin-glass behavior. 
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II. THE MODEL 



The model we study is a totally connected one with Hamiltonian 

' 



H = -l y £K ij S i S j , (1) 
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the spins Si (1 < i < N) are Ising spins taking the values ±1. Here we take the interaction matrix K to be a nearest 
neighbor and next nearest neighbor type interaction 

K = JJ T (2) 
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where the Jij are independent Gaussian random variables such that Jjj = and JijJki = SikSji/N. Here the over- 
line denotes the disorder averaging and we note that the matrix J here is not symmetric. The same model with J 
symmetric may be studied, however taking J non-symmetric considerably simplifies the analytical study of the model. 

In this paper we will consider two models. The positive temperature model where K — JJ T (we will refer to this 
as the K > model) and the negative temperature model where K = —JJ T (we will refer to this as the K < 
model) . Contrarily to what happens in the usual Sherrington-Kirkpatrick model, the model ([!]) is not invariant under 
the transformation K — > —K. Clearly if is a matrix which is positive definite and —K a matrix which is negative 
definite. 

For a symmetric matrix J taken from the Gaussian ensemble, the density of eigenvalues A is given by the Wigner 
semi-circle lawB 



i 
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a/4 - A 2 
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with A G [—2,2]. In corresponding the spherical spin modefl, at low temperature, the system minimizes its energy 
via a macroscopic condensation onto the eigenvector corresponding to A = 2. In the case of the SK model this 
condensation is not possible due to the discrete nature of the spins, however it seems reasonable to assume that the 
density p(X) for A ~ 2 plays an important role in the low temperature behavior. In our model, if J is taken to be 
symmetric, the density of eigenvalues K is clearly given by 



P±k(k) 
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where now, k € [0,4] in the +K or K > model and k £ [—4,0] for the —K or K < model. In fact one can 
show thato this is also the density of eigenvalues in the case of J non-symmetric. Hence in the positive case, the low 
energy region of the interaction matrix has the same form as the SK model and one would have in a spherical model 
a condensation onto the eigenvector of largest eigenvalue. In the negative case however the largest eigenvalues are 
at k = 0, but the density of eigenvalues now diverges in this region. In the corresponding spherical model it is easy 
to see that this eliminates the finite temperature phase transition. This happens in exactly the same way that the 
divergence of the density of occupation of the zero energy states for free bosons in two or less dimensions eliminates 
the finite temperature Bose Einstein transition. 

In the following treatment we consider the K > model although the same mathematical treatment we present 
here can be applied to the K < case on changing the sign of the inverse temperature (3. The Hamiltonian of the 
model can thus be written as 



H — — ^ ^ JikJjkSiSj 



(5) 



and the partition function for the model is therefore given by 



Z = Trst exp 



(6) 



In order to facilitate taking the disorder average we make a Hubbard Stratonovich transformation by introducing the 
auxiliary Gaussian spins x$ to obtain 



Z = Tr St , Xl exp 



P-TJSx- 

i'j 



where Tr x indicates the trace over the Gaussian spin x and is defined by 

rp [J , (3X\ 

Lr x = / ax \ — exp( ) 

/ V 2?r v 2 ; 



(7) 



(8) 



In the form (0) the model is that of an asymmetric Little model! but where one of the two sets of spins is Gaussian 
rather than Ising. Interest in the Little model aepse one decade ago, in the context of neural networks as the 
parallel dynamics of the standard Hopfield modellilJ coincides-.with the sequential dynamics in the Little modellill. 
The corresponding mean-field spin-glass model was studied ir£3 where it was shown that the equilibrium behavior is 
the same for both the SK and Little (for both asymmetric and symmetric couplings) models. Here we have verified 
numerically that, in our model also, a symmetric or asymmetric matrix leads to the same physical behavior. 
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III. THE POSITIVE TEMPERATURE K > MODEL 



We introduce replicas of both the Gaussian and Ising spins in order to average over the disorder via the replica 
method, obtaining 



Z n = Tr s *, x ° exp 



Trs a ,x a exp 



2N 



ij ab 



(3 2 N 



E QabPab 



where we have introduced the Ising and parabolic spin overlaps 

Qab = ^J2 S i S i and Pab=^J2 X i X i 



i z 

The trace over the Ising spins is accomplished using a delta function representation of the overlap constraint 



Tr s °5 



N 



Q ab - S i S i^j = (^p Tr S? / dA -" eX P ^ Y E A ^Qab - \Kb E S?S% 



where 



Z s = Tr sa exp E A ^S a S b ^j 



The same procedure is used for the Gaussian spins to yield 

extr S**[Q,P,A,T] 



In Z n 



N 



where S** [Q, P, A, T] is the saddle point action 



S** [Q, P, A, r] = L ^ QabPab + \ E A abQab + \ E V a b P a b 
ab ab 

+ Mz s ) + Mz x ) 



ab 



and 



Tr x aexp ^-lj2 T abX a x b ^j 



(9) 
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The saddle point equations with respect to P a b and Q a b allows the evaluation of the saddle values of the variables A Q 6 
and r afc 



dS** 

dPa~b 

dS** 

dQab 



= => Tab 



-p 2 Qab 



= => Kab = -l3 2 Pab 



leading to the reduced saddle point S*[Q, P] such that 



In [Z n 
N 



extr S*[Q,P] 



(16) 
(17) 

(18) 
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with 



S*[Q,P] = -^J2Q ab P ab + ln 



-Tr log (I -pQ) 



and where the Gaussian integral giving Z x has been evaluated. The saddle point value of P ab is given by 



thus leading to the new reduced action S[Q] such that 



In Z r 



N 



= extr S[Q] 



with 



\ (n -Tr(I- pQ)' 1 - Tr\n(I - PQ)) + lnTr s « exp f | £ [(/ - PQ)' 1 ]^ S a S b ^j 



S[Q] 



The saddle point equation from (123) is 



-^(I-(3Q)- 2 + ^(I-PQ': -(/- ^> -.1 = 



-1 ■ £ 

2 



where the elements of the matrix A are given by 

Tr s . S a S b exp (f £ ab [(/ - PQ)- 1 ] ab S a S' 



Tr S oexp(f E ab [(I-PQ)-\ b S a S b 



(19) 



(20) 
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We use the fact that the matrix (/ — PQ) should not be singular at the saddle point (otherwise the saddle point lies 
on a branch cut) to conclude 



Q = A 



(25) 



which is the physical saddle point equation for Q. We proceed by studying the replica symmetric from of this action 
with 



Q* b s = (1 - q)5 ab + qU ab 
where U ab = 1 for all (a, b). After some algebra one finds that 

s(q)= Um £i-LU = - In (1-/5(1 -q)) 

™— o n 2 2(1 — P(l — q)) A 



(26) 



iPq(-l + 2pq(l-q)) 



+ 



dz 
■J2i 



exp ( I In ( 2 cosh 



Pzy/q 



(l-P(l-q)) 



Taking the saddle point with respect to q above yields 
ds(q) P 2 (pq + p-1) 



dq 2 (1-/5(1-^)3 

There are two sets of possible solutions to Eq. (|2 



•?)) 



(27) 



(28) 



(29) 



and 
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* expf-^WW ^ ^f^ ) (30) 



2tt "V 27 V (1-/9(1-?)) 

However as mentioned above the solution (^) is unphysical. The solution (^) corresponds to the physical solution 
( p5| ) and only has the paramagnetic solution q = for T > T c = 2. Hence for T > T c = 2 one is (in the replica 
symmetric scheme) in a paramagnetic phase. 



The energy per spin _E obtained from (22) is 

1 



E = lim — 

n^O 2/3 



1 - -Tr(l - PQ) 



-1 



7? 



(31) 



where Q a b here is the solution to the saddle point equation. The replica symmetric energy per spin is thus 



— H RFf J (32) 

Hence for T > T c = 2, that is to say in the paramagnetic phase where one can be sure the the RS solution is correct, 
one has 

E = -W^W) (33) 

For T -C 1 a low temperature expansion of Eq. (BO) yields 



9 = 1- I -^ 7 = CM) 



With the replica symmetric ansatz, in the spin-glass phase, the ground state energy is 



1,2 , 2V2 



EBI=--(- + ^ + l) (35) 

As in the SK model if one computes the entropy corresponding to this replica symmetric solution within the spin- 
glass phase one finds a temperature where it becomes negative, indicating the need to break replica symmetry. The 
continuous nature of the transition however suggests that the underlying physics is the same as that of the SK model 
and the Landau expansion of S^Qaf,] has the same generic form as that of the SK model indicating a continuous RSB 
which sets in at T c . The introduction of NNN interactions thus does not qualitatively change the behavior of the 
model, in agreement with the previously discussed intuitive notion that (for a given spin type) the behavior of the 
density of the largest eigenvalues determines the nature of the transition. 

IV. THE NEGATIVE TEMPERATURE K < MODEL 

The Hamiltonian (Q) is not invariant under the transformation K — > —K. In this case the action for the replicated 
partition function is given by 

S[Q] = ^{n-Tr (I + PQ)- 1 -Tr ln(J + pQ)) + In Tr s « exp (™ f E + PQ)~ X ] ab ^^J ( 36 ) 

In the high temperature phase we expect the annealed approximation to be exact and find 

/ann = -^M2) + ^ln(l + /3) (37) 
The annealed entropy per spin S'ann is thus 



32 df 1,.,, , a , , P 



^ = -^t r H2)--Hl + P) + ^ TW) (38) 



The annealed energy is given by 
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1 

2(1 + /?) 



(39) 



The annealed solution can only be valid as long as Sann is positive. If there exists a Abound such that S ann (Abound) < 
for (3 > Abound, then Abound gives us an upper boundJbr the inverse temperature at which the equilibrium transition 
occurs, (3k- In the Adam-Gibbs-Di Marzio scenariaij Tk — 1/(3k corresponds to the temperature where the con- 
figurational entropy S'conf vanishes and the subindex K in (3 stands for the Kauzmann temperature which originates 
from the Kauzmann paradoxO which is based on the observation that an extrapolation of a high-temperature entropy 
cannot cross the low temperature solution. Sometimes Tk is also referred to as the ideal glass transition tempera- 
ture not to be confused with the experimental glass transition temperature T g in the context of finite dimensional 
glasses which turns out to be only a convention corresponding to an extremely large relaxation time as opposed to 
a diverging one. Because both entropies (S^nn and .S'conf ) are monotonically increasing functions of the temperature 
and S'conf < S anu this means that Abound > 0K is an upper bound of the ideal glass transition temperature Tk- One 
finds that Abound ~ 8.82 and hence T b ound = VA>ouud ~ 0.113. 

In the replica symmetric ansatz the equation for the Edwards- Anderson order parameter q is 



9 



dz 

7^ 



exp 



tanli ( 



(1 + /?(!-?)) 



(40) 



The examination of this equation shows that one can no longer have a continuous phase transition from q = 0. 
low temperatures one can show that a solution with q non zero exists and that 



For 



1 



T 



- - 1 

2 1 



0(T 2 ) 



(41) 



Numerically the inverse temperature (3^ at which q can become non zero is found to be (3? 5 « 29.3, thus giving 
j'RS 0.034. This is much lower than Tk above and hence the replica symmetric solution cannot eliminate the 
entropy crisis at Tbound- Clearly one must resort to a replica symmetry broken ansatz. Guided by the results of our 
RS calculation and the numerical simulations we make the REM (one step) like ansatz where the matrix Q is given 
by n/m matrices Q about the diagonal of Q of size m x m and is zero outside these blocks. The matrix Q takes the 
form 



Q = I(l-q)+qU (42) 
where / is the identity matrix and U is the matrix with each element equal to 1. We note the following results 

1 - !3qU 



{I + 0Q) 1 l + /?(l-g) /+ (l + /?(l-g)HJ-r HI - q ) + I ma l 



(43) 



Tr(l +/3QY 



m 
n 



Tr(l + (3Q)- 1 
1 



m — 1 



m \1 + /?(! - q) +/3mq l + 0(l-q) 



(44) 



Trln^l+ZJQ)- 1 ] = — Trln (1 + /3Q)" 1 



- (ln(l + (3{l -q) + (3qm) + (m - 1) ln(l + (3(1 - q))) 
m 



(45) 



and 



lnTr So expf| ^ [(I - (3Q)- 1 } ab S a S b j = £ lnTr s « exp ( t £ [(/ + pQ)~ x 

\ ab / \ ab 



n(3 



2(1 + 0(1 -q)) 



+ — In 

m 



dz 1 

—=exp(—-z 2 ) (2 cosh((3 za(q,m, (3)))' 
\y2n 2 



(46) 



G 



where 



a(q,m,8) = 

One therefore finds the one step action 

S[Q 1RSB ] 



V(l + P(l-q))(l + 0(1- q) +Pmq) 



(47) 



s(<7, m) 



lim 

1 

2 



-/3 2 g(l - q + mq) 



(1 + /3(1 -«))(! + )9(l-g)+m/3?) m 



ln(l+/3(l-g)+m/3g) 



- (1-1) In (l + ^(l-g)) 

TO 



H In (z(g, to)) 

to 



where 



One finds that 
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z(q,m) 



— s(q, m)= 8 a(q, to, 0) 



da(q, m, 0) 



1 



q 



dz 



dz 1 

J ex P<~2 z2 ) ( 2c °sh(/3za(g,™,/3))) r 



(1 — m) 
1 



(48) 



(49) 



z(q, to 

The saddle point of s(q,m) corresponding to the solution 

1 



exp(— — z 2 ) (2 cosh(0za(q, to, 0))) m tanh 2 (/3a(g, to, /3)z) 

27T 2 



is thus 



'7 = 



z(q, to) 



dz 1 

== exp(— -z 2 ) (2 cosh(/3za(q, m, /3))) m tanh (0a(q, m, /3)z) 



(50) 



(51) 



If one considers the case m = 1 this solution must give the same free-energy as the annealed free energy. However one 
can find a non zero value of q which signals a dynamical transition where the system becomes stuck in metastablc 
states of high free energy. Setting m = 1 in Eq. (O) gives 



exp 



±P 2 a(q,l,0) 2 



dz 1 

^ exp(— -z 2 ) cosh(/3za(g, 1, 0)) tanh 2 (/3a(q, 1, 0)z) 



(52) 



For small this equation has only the solution q = 0. However at 0d ~ 7.325 (Td = l/0d ~ 0.137) one finds a non 
zero value of q with q « 0.922. This transition corresponds to whatsis known as the mode-coupling transition in 
mode-coupling theories of the glass transition in their idealized versior£3. 

The precise way to locate both,transitions (Td and Tk) was suggested in a series of papers by Kirkpatrick, Thit.umalai 
and Wplynes and collaboratoral3 and later on applied to several models such as the random orthogonal modeO, Potts 
glassesE3 and mean-field quantum spin glassestj. The static transition is located by expanding the corresponding free 
energy f(q,m) about to = 1 and writing 

0f(q, to) = /3/ para + (1 - m)V(q) + O ((1 - to) 2 ) 

where / para = /(0, 1) is the paramagnetic free energy and 

df(q,m) ds(q,m) . 



(53) 



V(q) = -0- 



dm 



dm 



(54) 



The static transition Tk is given by the values of (0k where V(qu) = V'(qK ) = 0. The dynamic transition Td is given 
by the conditions V'(qd) = V" '(qd) — which correspond to the marginality condition and coincide with the solution 
found at to — 1 and reported above. 



The explicit form of the potential V is 

0q(l + 20) 



V(q) 



2(l + 8)(l + b(l-q)) 
1 



exp 



~-P'a(q,m,0)' 



+ -ln(l + /3(l-g)) - -ln(l + /3) 



dz 1 

exp(— -z 2 ) cosh(0za(q, 1, 0)) In (cosh(/3a(g, 1, 0)z)) 
2tt 2 



(55) 
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Using this to compute the static transition one finds that this gives Tk ~ 0.116 and qx ~ 0.985. The behavior of 
V(q) on lowering T is shown in Fig. ([!]). 

0.08 



0.06 



v (q) o.o4 



0.02 





0.2 0.4 0.6 0.8 1 

q 

FIG. 1. Behavior of the potential V(q) on lowering T. The different regimes are: T > Td (solid line), T — Td (dotted line), 
T K <T <T d (dashed line) and T = T K (dotted dashed line). 

Note that at Td the value of the potential V(qd) is finite and corresponds to the configurational entropy at 
that temperature. As T decreases below Td, the analytical continuation of the solution qd to the saddle point 
of Eq. (fig) increases while the value of rrid(T) decreases. The configurational entropy evaluated at T < Td is 
S C (T) = f3(f(qd(T), rrid(T)) — ipara(T)) and decreases while T decreases down to a temperature Tk where it vanishes 
and qK is the analytically continued value of (^(Tr-). 

Hence we see that in the case K < the physics of the problem is drastically altered by the interaction matrix, 
and that the model now exhibits the phenomenology of a (mean field) structural glass. 

V. NUMERICAL SIMULATIONS. 

We have verified the main predictions for the NNN-SK model with positive and negative temperature by doing 
some numerical simulations in both cases. For the positive temperature model we checked that the transition indeed 
occurs at T = 2. This has been done applying standard finite-size scaling techniques useful for investigating small-size 
systems. For the negative model we found that-strong freezing occurs at the mode-coupling temperature Td = 0.137 
as happens for other models such as the p-spinEa or the random orthogonal modelE-O. 

A. Some details of the simulations 

Simulations consist of standard Monte Carlo annealings using the Metropolis algorithm. The system is cooled 
down from high temperatures (typically twice the value of T c for the K > model and twice Td for the K < case). 
Annealing schedules are as follows: every AT = 0.2 for the K > model and every AT = 0.01 for the K < model 
the system is allowed to equilibrate over 1000 Monte Carlo steps (MCs) and statistics are collected during 10 5 MCs 
at each MCs. The sizes are small, N = 25, 50, 75, 100, but enough to locate the transition with some precision. The 
simulated range of temperatures are from T = 4 down to T = 0.2 for the K > model and from T = 0.3 down to 
T = 0.01 for the K < model. The number of samples where several thousands for all sizes. 
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Due to the long-range character of the interactions the dependence of the time needed to do a MCs grows quite 
fast with the size of the system (actually like iV 2 ). Therefore, for the statics, we had to limit our investigation to 
relatively small sizes. Moreover, a careful study of the relevant parameters for the transition (such as the kurtosis 
or the G parameter to be defined below) requires a large number of samples (this is especially true for parameters 
like G which measure sample-to-sample fluctuations). This last parameters the most successful example of what are 
referred to as order-parameter fluctuation parameters (OPF parameters)! 2 .]. 

Before showing the results let us mention that, while for the positive K model we achieved thermalisation in a 
range of temperatures in the vicinity of T c , for the negative K model thermalisation was hardly achieved due to the 
quite small acceptance rate for all the temperatures simulated. This behavior is due to the small value of the relevant 
temperatures of the K < model where the transition occurs (one order of magnitude smaller). Because the typical 
energy change for both models is the same, the Boltzmann factor is drastically reduced for changes which increase 
the energy in the K < model as compared to the K > case. This implies a very small acceptance rate for the 
K < models as compared to the K > case. In Fig. (|^) we show the acceptance as a function of T for the two 
cases. Note that the acceptance is nearly two orders of magnitude smaller in the negative model as compared to the 
positive model. Hence a good sampling of the configurational space for the negative model can be excluded. 





FIG. 2. Acceptance rate as a function of T in the relevant temperature ranges for the K > (right set of curves) and K < 
(left set of curves) models. For each set and from top to bottom we have N = 25, 50, 75, 100. Note that the acceptance rate is 
typically 10 times smaller in the K < model. 



B. The positive temperature K > model 



In this model there is a continuous RSB transition at T c = 2. The energy and the specific heat are shown as a 
function of T in Figs. (|J) and (Q). We also plot the result for the annealed expression only valid above the critical 
temperature. As we see the behavior of these quantities is similar to what is found for the SK model: the maximum 
of the specific heat occurs below T c . 



<) 





The transition can be well determined by looking at the B and G parameters as functions of the temperature. 
These parameters are defined by, 



2 



my 



G Am 2 -m 2 (57) 

<9 4 >-((<? 2 » 

In the infinite-size limit these parameters behave as: B{T) = B(T)9h(T c — T) while the behavior of G turns out 
to be simpler, G{T) — ^9h(T c — T) and hence transpires to be a better indicator for the transition. In Els. ( |5£j ) we 
show the kurtosis parameter B as a function of T for different sizes. In the same way as in the SK modelEd we find a 
crossing at a temperature close to T c = 2. This crossing turns out to be also present for the OPF parameter G 
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FIG. 5. Kurtosis B as a function of T for the model K > 0. From top to bottom we have N = 25, 50, 75, 100. The crossing 
temperature shifts with the size towards T c = 2. 




FIG. 6. OPF G as a function of T for the model K > 0. 
different curves cross close to T c = 2 where G c ~ 0.057. 



From top to bottom at low T we have N = 25, 50, 75, 100. The 



C. The negative temperature K < model 

As mentioned previously, we were not able to thermalise for this case close to the transition so we do not have good 
data for order parameters such as B or G. We show the results for the energy and specific heat in Figs. (0j8|) . 
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As we said before we did not succeed in thermalising the model at low temperatures to extract the behavior of the 
order parameters. Consequently we do not have good data for the kurtosis B (|56|) and the OPF G ((5?]). The only 
transition about which we may have some hints is the dynamical transition where dynamics becomes extremely slow. 
One could be tempted to interpret the crossing point for the specific heat at different sizes (Fig. (JsJ) ) as the signature 
of that dynamical transition. Despite of the fact that this interpretation seems reasonable we are not absolutely 
certain and so we prefer to leave this question open. 

More evidence for the dynamical transition can be obtained by doing annealing experiments for finite cooling rates 
and very large sizes. One of these cooling experiments is shown in Fig. (||) for a system of size N — 1000 was 
cooled down to very low temperatures staying for 10 5 MCs at each temperature and changing T by 0.01. The system 
gets trapped below a crossover temperature T* ~ 0.14 and starts to deviate from the high-T line. We plot two 
independent dynamical histories to demonstrate how the departure from the ergodic line depends on the dynamical 
history signaling the presence of non-ergodic effects in the dynamics (below T* the typical relaxation time has become 
much larger than the time the system remains at each temperature, i.e. 10 5 ). 
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FIG. 9. Cooling experiments in the model K < 0. The two vertical dashed lines limit the range of temperatures 0.13 — 0.14. 
The dotted line is the high-T result 

Note that in this dynamical study the static transition Tk remains completely hidden. At this temperature no 
divergence of a relaxation time can be observed since that time has already diverged above at the dynamical transition. 
The only way to show the existence of this transition is by estimating the configurational entropy and determining 
the temperature at which-that quantity vanishes. Such a task has been successfully undertaken in the case of the 
random-orthogonal modelE3 which displays a very similar behavior to the K < model and requires the use of 
projection techniques, in the spiritj-af those introduced by Stillinger and Weber, for the study of the potential energy 
landscape in the context of glasses^!. 

VI. CONCLUSIONS 

In this paper we introduced a solvable model of a spin-glass corresponding to the mean-field version of the nearest 
and next nearest neighbor spin-glass model. The main interest of the model is that it shows how correlations in the 
couplings may completely change the character of the spin-glass transition from continuous replica symmetry breaking 
to a one-step transition. We analyzed two cases, the positive model K > with a phase transition and behavior similar 
to that of the SK model and the K < model with ideal mode-coupling behavior and a phase transition similar to the 
ROM or to the p-spin model. The different character of the transition in both models can be ascribed to the different 
behavior of the eigenvalue density distribution close to the maximum eigenvalue or threshold. For the K > the 
distribution vanishes at the threshold (i.e. A + = 4) and the equilibrium distribution condensates around configurations 
in the vicinity of the maximum eigenvector. For the K < the eigenvalue density diverges at the threshold (i.e. 
A_ =0). In that case, the equilibrium configuration condensates around an extremely large number of eigenvectors 
with eigenvalue close to 0. Hence, the phase space splits into a very large number (exponentially large with the size 
of the system) of ergodic components or phases, these phases corresponding to different eigenvectors with eigenvalue 
very close to the threshold which extensively contribute to the configurational entropy. Obviously, eigenvectors of 
the coupling matrix never coincide with possible configurations of the Ising system, so this argumentation must be 
taken only as a rough picture. Nevertheless, the idea that the type of eigenvalue distribution determines the character 
of the transition seems quite intuitive. Actually, if one considers spherical instead of Ising spins then the transition 
disappears for the K < model but persists for the K > case. For the K < model the transition disappears 
because there is no longer a vanishing of the configurational entropy (the classical nature of spins allows for a negative 
entropy). The same mechanism occurs in the ROM where the one-step transition with Ising spins disappears in case 
of spherical spins. In this case, the density of eigenvalues is given by p(X) = |<5(A — 1) + hS(X + 1) showing that the 
density diverges in the vicinity of the maximum eigenvalue. 

What is the real connection between the density of eignevalues and the order of the transition? The mechanism 
which makes a continuous transition become first order (in the order parameter) in spin-glass theory is pMnlL known 
and based on the role of the so called complexity or (sometimes misleadingly) configurational entropj£3't3t3. In 
mean- field theory the complexity S C (F,T) at a given free energy and temperature corresponds to the logarithm of 
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the number of TAP free energy minima solutions with that free energy and temperature. The complexity defines the 
so called potential <f>(F,T) = F — TS C (F,T). This quantity display two types of different behaviors at low enough 
temperatures. For models with continuous replica symmetry breaking $ has a single minimum at the equilibrium free 
energy value F cq (T) while for models with a one-step scenario this function may display a minimum at a threshold 
value F*(T) higher that the equilibrium free energy F cq (T). The difference 0(F*(T) — F cq (T)) defines the complexity 
at that temperature S*(T) — S C (F* (T),T) which is positive for one-step models but vanishes in continuous models. 
Because at high free energies $>(F, T) ~ F it is clear that the abundance of metastable solutions makes S c always larger 
raising the possibility that the potential $(F, T) displays a minimum above the equilibrium free energy. Therefore, a 
divergence of the eigenvalue density at the lowest value goes in the appropriate direction of generating a minimum in 
$ and hence making the transition become first-order in the order parameter. 

As for further directions of interest we mention the investigation of the existence of a model interpolating between 
the K > and K < models. This would be very interesting to understand better the mechanism which makes the 
spin-glass transition change from a continuous one to a discontinuous one, as well as to connect this change with the 
behavior of the density of eigenvalues close to the threshold. In this direction it would be also interesting to investigate 
in a general way the type of phase transition for a generic eigenvalue distribution using analytical techniques such as 
those developed for the ROM modeO. 

It would also be interesting to investigate the behavior of the different terms in the TAP expansion to see how 
this change of behavior occurs and how this is related to the geometrical properties of the free energy landscape. 
Note that the information contained in the density of eigenvalues is closely related to the topological features of the 
energy landscape. Therefore, it is reasonable to think that the nature of the transition in these models arises from 
the geometrical properties of the energy (and hence, the free energy) landscape. 

Finally, we would like to investigate the behavior of this model on a finite dimensional lattice. In the simplest way 
this can be done taking J to be nearest- neighbor Gaussian matrix and defining K according to (|^) . It would be very 
interesting to see how the two universality classes describing the mean-field behavior are modified in the presence of 
activated processes and how they manifest themselves in the dynamical properties of realistic systems. This could 
give additional understanding about why some models behave in one way or another and how topological properties 
of the landscape influence the dynamical behavior of real systems. 
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